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Abstract . Since the force of gravity is inverse- 
ly proportional to the square of the distance be- 
tween two masses, it is obvious that the effect of 
small gravitational anomalies on the motion of an 
orbiting satellite will increase as the satellite 
altitude decreases. However, at very low altitudes, 
the effect of atmospheric drag results in drastical- 
ly reduced orbit lifetimes and considerable uncer- 
tainty in satellite motions. The concept suggested 
herein employs a Dis turbance compensation System 
(DISCOS) on each of a pair of satellites at very 
low altitudes to provide refined measurements of 
the earth's gravitational field. The DISCOS main- 
tains the satellites in orbit and essentially elim- 
inates motion uncertainties due mostly to drag and 
to a lesser extent from solar radiation pressure. 

By a closed-loop measurement of the relative range- 
rate between the two low satellites, one can de- 
termine the earth's gravitational field with a 
considerably greater accuracy than could be ob- 
tained by tracking a single satellite. 

Introduction 

Since the advent of artificial satellites of the 
earth, the science of geodesy has advanced with re- 
markable rapidity. Without employing orbiting sat- 
ellites, our knowledge of the earth's gravity field 
could not be determined to nearly the level to 
which it is now known. This improved knowledge of 
geodesy is accomplished essentially by studying the 
long period motions of orbiting satellites, partic- 
ularly satellites at quite moderate altitudes. The 
lower the altitude of the satellite, the more pro- 
nounced the gravitational effect of the earth's 
gravity field, particularly the higher order har- 
monics. Although it is obviously desirable to go 
to still lower altitudes, the uncertainty of the 
along track force caused by air drag results in a 
severe degradation of the accuracy to which the 
satellite motion can be measured. 

To determine the earth's gravity field, with 
greater precision it is necessary to make increas- 
ingly more accurate measurements of the satellite's 
motion. A basic limitation in the use of a single 
satellite to determine the earth's gravity field is 
that it is very difficult to measure the minute 
difference in velocity of the low orbiting satel- 
lite due to small gravitational effects in the 
presence of the overwhelming 7 kms/sec of satellite 
orbit velocity. 

Suggested herein is the concept of providing a 
drag-free capability to each one of a pair of sat- 
ellites at very low altitudes. The principal ad- 
vantages of such a system are that l) A drag free 
system allows the satellite to go to extremely low 
altitudes where they are most sensitive to gravi- 
tational anomalies and meaningful velocity measure- 
ments can be accomplished because the satellite is 

Proc. of the 9th GEOP Conference. An International Symposium on the Applications of 
Geodesy to Geodynamics, October 2-5. 1978. Dept, of Geodetic Science Rept. No. 280. The 
Ohio State Univ., Columbus. Ohio 43210. 


(ONLY ONE AXIS IS SHOWN) 



Fig. 1 Illustration of DISCOS as used on the Triad Satellite 


free of the effects of atmospheric drag; 2) Pre- 
cise measurements of the relative velocities of 
two satellites in identical orbits is vastly sim- 
pler than the measurement of the absolute velocity 
of a single satellite. This is because it is less 
difficult to directly measure a small quantity 
rather than to make a measurement of a small quant - 
tity in the presence of a very large quantity. 

DISCOS Concept and Performance 

A Dis turbance Compensation System (DISCOS) has 
been shown to perform satisfactorily in orbit on 
the TRIAD satellite. An illustration of a single 
axis of the three-axis DISCOS system that was used 
on TRIAD is shown in Figure 1. The TRIAD satellite 
was extended in the vertical direction to achieve 
gravity-gradient stabilization. As such, it pre- 
sented a considerable cross-sectional area on which 
the atmosphere acted to create drag. At the center 
of the satellite was a dense spherical mass called 
the proof mass which was contained within a spher- 
ical cavity that shielded the proof mass from the 
effect of external disturbance forces such as solar 
radiation pressure and air drag. Contained within 
the spherical cavity at the center of the satellite 
was an imaginary spherical cavity known as the dead 
zone. When the effect of drag was sufficient to 
push the satellite backward until the dead zone 
touched the surface of the proof mass, the aft 
thruster fired thereby causing the satellite to 
move forward relative to the proof mass. By this 
means, the proof mass was allowed to fly in an or- 
bit determined purely by gravitational forces, and 
the satellite, because of the DISCOS system, was 
constrained to follow the motion of the proof mass. 
In Figure 2 is shown a particular 12 minute period 
of data collected from the TRIAD satellite. In 
this figure is shown the along-track position on 
the proof mass expressed in millimeters as a func- 
tion of time from the start of a pass of the satel- 
lite. During the time of this pass, the principal 
disturbance force on the satellite was drag, so that 
at approximately one minute of time, the integrated 
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drag force was sufficient to cause the proof mass to altimeter is shown which might provide improved 
touch the imaginary dead zone. This is a consequence data on ocean surface topography which of course 
of 1.0 mm of motion of the proof mass from the geo- could contribute materially to our knowledge of 


metric center of the cavity. At that point, which 
is shown in Figure 2 as the first of three "Aft 
Thruster Fires" indications, the impulse of the aft 
thruster causes a velocity of the satellite rela- 
tive to the proof mass such that proof mass has an 
apparent motion back towards the geometric center 
of the cavity. However, due to the fact that the 
drag forces are continuing at this time, the proof 
mass never reaches the center of the cavity. Four 
minutes after the first of the thruster firings, 
the proof mass again touches the dead zone. The 
parabolic shape of the displacement curve of the 
proof mass as a function of time is an indication 
of an essentially constant drag force during this 
four minute period. From the times of 5 minutes to 
approximately 11 minutes, the satellite is again 
pushed away from the proof mass but this time the 
different shape of the essentially parabolic curve 
indicates that the drag force at that time was 
smaller. Finally we see at 11 minutes the third of 
the firings of the aft thruster. This data is typ- 
ical of what occurs in a DISCOS system as the fir- 
ing of the thrusters are used to keep the proof 
mass centered in the satellite's spherical cavity. 

Application of DISCOS to Satellites 
at Very Low Altitude 

The TRIAD satellite containing the first three- 
axis DISCOS flew at an altitude of approximately 
700 km and used three pounds of cold gas (Freon) 
propellant to achieve 18 months of drag free oper- 
ation. If one is to maintain a spacecraft at a 
much lower altitude or if one is to maintain drag 
free operation for a longer period of time more 
total impulse from the propellant is required. 

This can be accomplished either by providing a 
greater quantity of propellant and/or by providing 
a much higher specific impulse as compared with 
cold gas. In Figure 3 is shown the concept of a 
satellite which could be one of a pair used in a 
drag free lo-lo system for obtaining more refined 
measurements of the earth's gravity field. The 
satellite shown in Figure 3 has a DISCOS system at 
its center mass and two very large tanks contain- 
ing propellant (e.g., hydrazine) located symmet- 
rically on each side. It is estimated that for a 
satellite of this size, 2000 kilograms of hydra- 
zine could be readily provided. Figure 3 also 
shows the concept of simple angular momentum fly- 
wheel that could be used to stabilize the satellite 
in roll while aerodynamic forces could be used to 
stabilize the satellite in pitch and yaw. A radar 
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the earth's gravity field and ocean dynamics. 
Furthermore, this concept shows a vector and sca- 
lar magnetometers which are indicative of another 
class of measurement which is best accomplished 
with satellites at very low altitudes. 

In Figure 4 is shown the propellant usage for 
such a satellite in the altitude region between 
100 and 250 kms. The propellant usage in Figure 
!(■ is expressed in kilograms /month for a satellite 
using hydrazine as the DISCOS system propellant 
and having a frontal area of 1 square meter. Two 
curves are shown in Figure 1*; one for a solar min- 
imum and the other for a solar maximum. If a sat- 
ellite such as that shown in Figure 3 can contain 
2000 kilograms of hydrazine, it will provide ap- 
proximately 2 months of drag-free operation at a 
125 km and approximately 20 months operation (at 
a solar minimum) at an altitude of I 5 O km. 

A consideration which makes the satellite con- 
cept as shown in Figure 3 reasonable from a cost 
standpoint is the advent of the Shuttle as a launch 
vehicle. The Shuttle is particularly well suited 
to launch satellites that are extremely heavy but 
that require very low altitude for their operation. 
If longer times in orbit are required, the Shuttle 
astronauts could rendezvous with the pair of drag- 
free satellites and refuel these spacecraft in- 
orbit . Although there would be times when it is 
desirable to fly the satellite at altitudes as low 
as 125 to 150 kilometers, most of the time a con- 
siderable improvement in the knowledge of the 
earth’ 8 gravity field would result from flying the 
satellites at altitudes in the region of 250 kilo- 
meters. At such an altitude, even at the solar 
maximum, 2000 kilograms of hydrazine would provide 
approximately 200 months of in-orbit operation 
without refueling. This mission duration time 
would be satisfactory for significantly improving 
the accuracy to which we know of the earth's 
gravitational field. 

Another effect that must be considered in al- 
lowing spacecraft to fly at extremely low 'altitudes 
is the aerodynamic heating that results from fric- 
tion of the high speed spacecraft with the earth's 
atmosphere. In Figure 5 is shown the temperature 
at the stagnation point which would be obtained at 
the leading surface of a satellite that was flying 
at altitudes between 100 and 250 kilometers. The 
assumptions here are that the front surface has an 
infra-red emmisivity of 0.8 and there is no energy 
input from the sun. From these curves (which show 
solar maximum and solar minimum), it can be seen 
that at an altitude of 125 kilometers, the skin 
temperature at the stagnation point under the con- 
ditions described above is below 300^* C. Such a 
temperature could readily be withstood even by 
aluminum. When one rises to an altitude of 150 
kilometers. Figure 5 shows that the temperature 
of the front surface is not very much different 
from room temperature and therefore would not 
compromise the design of a spacecraft. In essence, 
Figure 5 shows us that there is no significant 
temperature barrier to operating a single satellite 
or a pair of satellites at extremely low altitudes; 
i.e., above 125 tan. 





Fig. 4 Propellant usage for satellites at very low altitude. 


Analysis of Lo-Lo System of Satellites 

In this section, an approximate expression is 
developed for the relationship between the gravita- 
tional environment and the relative range-rate be- 
tween the satellites of a lo-lo system. Through 
judicious approximations an analytic formulation 
can be obtained. This is useful to enhance concep- 
tual understanding and to provide a basis for eval- 
uating and interpreting the more detailed numerical 
simulation which follow. 

Assume that the two satellites In Figure 6 are 
in identical circular orbits but separated in phase 
by the angle a and have mean anomalies where i = 
1 or 2 which identifies each one of the two space- 
craft. Let T be the kinetic energy per unit mass 
and U the gravitational energy per unit mass. If 
the rotation of the gravitational field is neglected 
it is possible to introduce the principle of conser- 
vation of energy 



T + U = K (1) 

where K is a constant. Let L^, H^, be the spa- 
tial perturbations of the spacecraft from the cir- 


Fig. 6 Geometry of lo-lo system of satellites and 
definition of Lj, Cj, Hj coordinate systems, 
i = 1 to 2. 
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cular reference orbit consistent with the reference 
systems in the figure. Then 


where F is soley a function of the inclination 


T = A 
1 a 


(v, + if + f + f 


( 2 ) 


where v is the nominal circular velocity. Let 
o 


U - V - 


(3) 




v L. + v - - + di . 2 + H . 2 + C 2 ) 
ox ir l i i / 


= K 
00 


The bracketed term is of second order since L^, 
Hi, Ci « v Q and can be neglected. 

For closed circular motion 




(5) 


and 


V o = (a) (6 > 

Substituting Eqs. ( 5 ) and (6) into Eq. (k) give 


- -V 1 V(Mi) . 


(7) 


and 


jjmp 


imp 


lm 


-S. 


where V is the perturbing potential per unit mass, 
k is the gravitational constant and a the radius to 
the satellite. Substituting Eqs. (2) and ( 3 ), Eq. 
( 1 ) gives 


lm 


2 “ 


lm 


(l~ m)even 

(l -m)odd 
(i*m)even 

(i-m)odd 


Cos 


Sin 


(^-2p)(MF u) ) + m(n- 0 ) 

(13) 

(j6-2p)(m+uj> m(n-e) 


where R is a scaling distance, u, the argument of 
perigee, f) the right ascension of the ascending 
note, g the Greenwich sidereal time and l and m 
the order and degree respectively. For a polar 
satellite (q-q) is constant since the rotation of 
the earth has been neglected here. Consequently 
Vj is a harmonic function with arguments (i- 2p) 
(M+ca). 

Partitioning the potential in components 
permits Eq. (10) to be rewritten as 


where 


P. = -v 
2 o 


-1 


p = £ p„ 
2=2 * 


VV - 


(It) 


(15) 


The relative range-rate between the two spacecraft 
in the same plane is 

p = (^ + H 2 ) sin | + (^ - L g ) cos | . ( 8 ) 

For small separation distances, i.e., a small, it 
is appropriate the range-rate be 

P ” ^2 (9) 


For the lo-lo configuration 

Mj_ = M + | , Mg = M - | (16) 

so that 

WV - S imp (M 2 ) " 2 sin^f^i^pCM^^e) 

(17) 


which using Eq. (7) becomes 

p = -v 0 ‘ 1 (V.fttj.) - v(m 2 )). (10) 

Under the assumptions stated, this equation demon- 
strates that the relative range-rate is a measure 
of the difference in the perturbing potential. This 
result has been previously derived by Wolff [1969] 
and Comfort [1973] • Now, let the perturbing poten- 
tial be expressed by 


V = £ V. . (11) 

2=2 i 

For a near circular orbit Kaula [1966] gives 


where S^ m p is the derivative of S^ m p with respect 
to its argument and u>=un=u>2> fl = On = Og * Substi- 
tuting Eqs. (12) and (17) into Eq. (15) gives 



P 2 

£ £ 


m=o p=o 


2 sin 


F s 

2 imp imp 


(18) 


The bracketed term represents a function dependent 
on the separation of the two spacecraft that modu- 
lates the amplitude of the range-rate. This term 
demonstrates the tuning capability of the lo-lo 
configuration. For a particular value of i say i 0 , 
Eq. (18) shows that frequencies i 0 , i c -2, i 0 -k, • • • 
times (tttu>) exist where contributions to the range- 
rate with frequency 2 0 (M+uj) results from coeffi- 
cients Cj m and S^ m such that 



2 

£ 

p=o 


imp 


S imp (uj ' M > n » 9) 


( 12 ) 


2 = i Q + 2 p p = 0 , 1 , 2 , •• • (19) 
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In this context, Eq. (19) shows that by selecting 
a separation 

a = nlgstil n = 0,1, 2, • • • (20) 

•®o 


so that Eq. (27) can be rewritten as 



Substituting this into Eq. ( 15 ) gives 


the contributions to the range-rate from coeffi- 
cients C im and S, m 

4 = i Q + 2p p = 0,1,2, ••• (21) 


-1 


_R_ 

1 - 1 ' 


4+1 


Ag, (Mj+u>) - Ag, (M 2 +uo) 


(31) 


will be doubled. Conversely, all coefficients 
such that 


as the expression for range-rate. The maximum 
amplitude for will occur at the spacing given 
by Eq. (20) where 


4 = 2mU 0 + 2p) n = 0,1,2, •• • (22) 


will contribute nothing to the signal. This tuning 
capability is illustrated in Figure 7- Another way 
to conceptualize tuning is to consider the separa- 
tion to be such that the variations in the gravita- 
tional field are either negatively or positively 
correlated. Enhancement in the signal, will occur 
with the former and destruction with the latter. 

An estimate of the amplitude of the range-rate 
for a specific gravity anomaly can be obtained as 
follows. The gravity anomaly Ag at the spacecraft 
can be written in terms of the disturbing potential 
V as 


Ag 


21 

dr 


r=a 



(23) 


Substituting Eqs. (ll) and (12) for V gives 


Ag - 


R \* +2 L * 


R 4=2 
Let Ag, be defined by 


ECi- 1 ) ft £ £ 


m=o p=o 


^mp imp 


Ag = £ Ag 
4=2 1 


where 


Ag 4 R 2\a 


,4+2 4 4 

E £ F, S„ 

„ 4mp 4mp 
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Comparing Eqs. (12) and (26), gives 


(24) 


(25) 


( 26 ) 
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k 4+l 
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+ -). 


(32) 
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Fig. 7 Effect of separation distance on velocity amplification. 


This equation has been used as the basis to deter- 
mine the accuracy in the range-rate measurement 
that is necessary to recover a specified gravity 
anomaly. Results are presented in Figure 8 where 
the curves represent the amplitude and wavelength 
of the gravity anomalies that produce a one way 
range-rate of the specified instrument accuracy. 

At the longer wavelengths a small angle approxi- 
mation for optimum separation becomes less valid. 

A recent study by Goldfinger [1978] have shown 
that for harmonics of order two even at the opti- 
mal spacing the error is less than a factor of 2. 

Corruption of the signal by errors in the 
knowledge of the initial conditions has been dis- 
cussed by Fisacane [1978] . There, a computer sim- 
ulation study show that the effects of the initial 
condition errors can be effectively eliminated by 
high pass filtering. 


Since the interest is in recovering a gravity 
anomaly whose magnitude is specified at the sur- 
face, let 

A^ = Ag, 

From Eq. (26), it follows that 

/ r \4+2 — 

Ag £ = (-) Ag, (29) 


Velocity Measurement Concept 

The relative velocity measuring instrumentation 
concept is indicated in Figure 9. The basic method 
can be described by thinking of the second satel- 
lite as a pure transponder whereby the tone (Nfj_) 
transmitted by the first satellite is simply sent 
back with a doppler offset resulting from the rel- 
ative motion between the satellites. The same 
doppler effect would also be experienced on the 
return trip resulting in a received frequency of 
Nf^(l - p/c)^ at the first satellite; where p 


( 28 ) 

a=R 
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Fig. 8 Anticipated accuracy for recovery of the gravity field 
based on relative velocity measurement accuracy of 
10“* mm/sec. 
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Fig. 9 System for determining relative velocity of a 
lo-lo satellite pair. 


FREQUENCY 
ANTENNA BEAM WIDTH 
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TRANSMITTED POWER 
RANGE 

SIGNAL-TO-NOISE RATIO 
INTEGRATION TIME 
la PHASE JITTER 
la VELOCITY NOISE 

Fig. 10 An example of lo-lo 


= 100 GHz 

= T x 10” 

= 18mm x 90mm (0.7" x 3.5") 
= 10 dB 

= 1.0 Hz 

= 2.0 mW 

= 112 km 

= 45 dB 

= 5.0 SECONDS 

= 0.06 DEGREES 

= 10 -4 mm/sec 

closed-loop link parameters. 


represents an average velocity over the measure- 
ment interval. If the return signal is then 
differenced with the original transmission freq- 
uency, the velocity component, approximately 
2Nf^ p/c, will be separated for measurement. The 
actual implementation incorporates a frequency 
translation in the second satellite to avoid 
oscillation difficulties. However, the implemen- 
tation provides for removal of both oscillator 
offset errors thereby providing the same final 
result. 

If both oscillators were properly characterized 
as having constant but unknown frequencies, i.e., 
having only bias errors, the technique just dis- 
cussed would fully remove all oscillator errors. 
When time varying errors are included, their ef- 
fect can in concept be removed by matching the 
delays into the differencing circuits. The delay 
shown in the figure between the f]_ oscillator and 
the return mixer is intended to carry out this 
function. Removal of time varying effects of the 
f 2 oscillator requires that the two delays between 
the in oscillator and the summing circuit in the 
second satellite be matched. Assuming that com- 
plete oscillator errors are removed by the imple- 
mentation, the remaining error sources are due to 
phase variations in uncommon circuits, propagation 
effects and thermal noise limitations in the link. 

Assuming all other effects can be made negligi- 
ble, the signal-to-noise performance shown in Fig- 
ure 10, can be realized with quite reasonable an- 
tenna dimensions and transmission powers. While 
the system is not fundamentally limited to this 
precision, a high degree of circuit phase stability 
will be required to achieve the indicated perform- 
ance. Prior to further evaluation of the instru- 
mentation limitations, the 10 _1 * mm/sec noise value 
is being considered as the practical measurement 
limit . 


Conclusions 

The advent of the shuttle launch system makes 
it possible to launch a pair of satellite contain- 
ing several thousand kilograms of propellant into 
a low altitude orbit. 

These large quantities of propellant can be used 
in a DISCOS system to maintain a satellite in orbit 
and to free the satellite from the otherwise over- 
whelming disturbance effects of atmospheric drag. 

A pair of such satellites at very low altitudes 
can greatly refine our knowledge of the earth's 
gravitational field. 

Acknowledgments . The authors want to thank 
Thomas Thompson and Andrew Goldfinger of the 
Applied Physics Laboratory for their assistance 
in helping in the preparation of this report. 

References 

The Staff of the Space Development Department, 

The Johns Hopkins University Applied Physics 
Laboratory, and The Staff of the Guidance and 
Control Laboratory, Stanford University, "A 
Satellite Freed of All but Gravitational Forces: 
TRIAD I, J. of Spacecraft and Rockets. Vol, 88, 
No. 9, pp. 637-644, 1974. 


218 





Wolff, M., "Direct measurements of the Earth's 

gravitational potential using a satellite pair," 
Journal of Geophysical Research, vol. 74. no. 

22, 15 Oct. 1969. 

Comfort, G., "Direct mapping of gravity anomalies 
by using Doppler Tracking between a satellite 
pair," Journal of Geophysical Research, vol. 78. 
no. 29, 10 Oct. 1973. 

Kaula, W. M., Theory of Satellite Geodesy. Blaidsdell, 
Waltham, Mass., 1 966. 

Goldfinger, A., Private communication, November 1978. 

Pisacane, V.L., "Filtering of Satellite-to-Satellite 
Observations: Low-Low Configuration," presented 

at the American Geophysical Union 1977 Spring 
Annual Meeting, Washington, D.C., May 30 -June 3> 
1977; (Acta Astronautica, in press). 


219 



